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Abstract. We present WeCAPP, a long term monitoring project searching for microlensing events towards M31. 
Since 1997 the bulge of M31 was monitored in two different wavebands with the Wendelstein 0.8 m telescope. In 
1999 we extended our observations to the Calar Alto 1.23 m telescope. Observing simultaneously at these two sites 
we obtained a time coverage of 53 % during the observability of M31. To check tho usands of frames for vari ability 
of unresolved sources, we used the optimal image subtraction method (OIS) bv lAlard fc LuptonI il998h . This 
enabled us to minimize the residuals in the difference image analysis (DIA) and to detect variable sources with 
amplitudes at the photon noise level. Thus we can detect microlensing events with corresponding amplifications 
y4 > 10 of red clump giants with A/i — 0. 
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1. Introduction 

In the last decade microlensing studies proved to be a 
powerful tool for searching baryonic dark matter in the 
Galactic halo. 

Several grou p s lik e the MACHO colla b oratio n 
llAlcock et all Il993^■ OGLE (Udalsk i et all Il993l) 
ERO S llAiiboiirg et all Il993l) and DUO llAlard et all 
Il99,'^^ followed the suggestion of lPaczvnskil l|l986l) and sur- 
veyed millions of stars in the Large and Small Magellanic 
Clouds (LMC, SMC) and in the Galactic bulge for vari- 
ability induced by gravitational microlensing. Although 
all of them discovered events compatible with gravita- 
tional lensing by MACHOs 



distribution in the Galact ic halo (e.g. iLasserre et all 
2000t lEyans fc Keriii j. |2000l and referen ces therein). 
Crottj ()1992|) and iBaillon et all l)l993(l suggested to in- 
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elude M 31 in future lensing surveys and pointed out that 
it should be an ideal target for these kind of experiments. 
In contrast to microlensing studies towards the LMC and 
the SMC, which are restricted to similar lines of sight 
through the Galactic halo, one can study many different 
lines of sight to M31, which allow to separate between 
self- lensing and true MACHO events. 
Since the optical depth for Galactic MACHOs is much 
greater towards M31 than towards the LMC, SMC or 
the Galactic bulge one expects event rates greater than 
in previous lensing studies. Furthermore M 31 contributes 
an additional MACHO population as it possesses a dark 
halo of its own. Thus, three populations may contribute 
to the optical depth along the line of sight: MACHOs 
in the Galactic halo, MACHOs in the halo of M31 and 
finally stars in the bulge and the disk of M31 itself, a 
contribution dubbed self-lensing. 

The high inclinat ion of M31 (77°) 
l|Walterbos fc Kennicutl . Il987(l produces a near-far 
asymmetry of the event rates. The near side of the M31 
disk wi ll show less events than the more distant one 
(ICrottsL Il992l) . Since Galactic halo-lensing as well as 
self-lensing events will not show this feature, a detected 
asymmetry will be an unambiguous proof for the existence 
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of M31 MACHOs. 

As most of the sources for possible lensing events 
are not resolved a t M3 1's distance of 770 kpc 
I Freedinan fc Madorel Il990() the name 'pixellensing' 
l GoiildLll99ir was adopted for these kind of microlensing 
studies. In the mid nineties two projects started pixel- 
lensin g surveys towards M3 1, AGAPE (lAngaj jeLali 
Il997l) and Columbia/VATT llTomanev fc CrottsL jwm 
First candidate events we re reported ijAnsari et all Il999t 
ICrotts fc Toma,nev[ll99fil) but could not yet be confirmed 
as MACHOs. This was partly due to an insufficient time 
coverage which did not rule out variable stars as possible 
sources. 

1999 two new pixellensing projects, P O INT-A GAPE 
i Kerins fc the Point- Agaoe CoUaborationl l2000(l . who 
reported recently a first candidate microlensi ng event 
l)Auriere et all l200lj) . and MEGA l|Crotts et a l.. 1999), 
the successor of Columbia/VATT, began their system- 
atical ob servations of M31. Another pro ject, SLOTT- 
AGAPE lISlott-Agane CoUaborationl 1 19991) will join them 
this year. 

The Wendelstein Calar Alto Pixellensing Project started 
with a test and preparation phase on Wendelstein as 
WePP in autumn 1997 before it graduated after two 
campaigns to WeCAPP in summer 1999 by using two 
sites for the survey. 

In this paper we will give an introduction to the project 
including information about the data obtained in three 
years and our reduction pipeline. In Sect. 2 we briefly 
discuss the basic principles of pixellensing. In Sect. 3 we 
will give an overview of the project including information 
about the sites used and the data obtained during 
WeCAPP. Section 4 refers to our data reduction pipeline 
and describes how light curves are extracted. In Sect. 5 
we show first light curves and Sect. 6 summarizes the 
paper. 

2. Pixellensing 

In microlensing surveys of uncrowded fields a resolved star 
is amplified by a function A which can be measured from 
the light curve and which yields direct information of the 
lensing parameters (Paczynski, 1986): 

o,2 I o 



u\/ v? + 4 



(1) 



with = l/^E, the inverse Einstein ring crossing time, 
uo the angular impact parameter in units of the angular 
Einstein ring radius and tg the time of maximum am- 
plification. 

The angular Einstein ring radius ^e is connected with the 
physical Einstein ring radius and the properties of the 
gravitational lens by 



rE _ / 4GM A A. 



(2) 



with the mass of the lens Af , and the distances between 
the observer and lens D\, observer and source -Ds, and 



lens and source As- According to Eq. ^a lensing event 
will produce a symmetric and, as lensing does not depend 
on the colour of the source, achromatic light curve (see 
also Sect. 3). 

In more crowded fields it is not possible to determine 
the amplification A unambiguously because many un- 
resolved sources may lie inside the solid angle of the 
point spread function Opsf of a bright source. If one of 
these unresolved sources is amplified, the event could er- 
roneously be attributed to the brigh t star, which r esults 
in a strong amplification bias liHanl 199 ^AlardL Il997l: 
IWozniak fc Paczvnskl Il997t iColdbergl Il99?j) . 
Similarly, in pixellensing studies the true amplification 
cannot be determined because many stars fall in one res- 
olution element. 

Nevertheless one can still construct a light curve Ai^(<) 
by subtracting a reference frame from the image, in which 
a lensing event takes place: 

b.F{t) = F(t) -B = F,- \A{t) -\\-B = F,+ B,,, (3) 

with Fi the flux of the source after or before lensing and 
Aes = ^j^i being the residual flux from unlensed 
sources inside fipsF- As there are many sources contribut- 
ing to the flux inside rJpsF, pixellensing events will only 
be detected if the amplification A or the intrinsic flux A 
of the lensed source are high. However, this difficulty in 
identifying events is balanced by the large number of pos- 
sible sources for le nsing ev ents in the bulge of M31. 
As pointed out bv ICouldl lll996.) . the main difference be- 
tween classical microlensing and pixellensing consists in 
the fact, that in the latter case the noise within flpsF is 
dominated by unlensed sources and therefore stays vir- 
tually constant during an event. As only events with a 
high amplification can be detected, the Einstein timescale 
^E is not a general observable in pixellensing. Therefore 
^FWHM which describes the width of a lensing light curve 
at half of its maximum value is the only timescale one 
is able to measu re with a. certai n accuracy. Based on a 
previous work of lGondolrJ ^2^, who pointed out that 
the optical dep th towa rds M31 can be estimated with- 
out knowing fo,. lBaltz~fc Silk ( 200Q1 showed how the mea- 
surement of a moment of the light curve ta-n permits to 
calculate the Einstein time t-E for a particular lensing 
event. Furthermore principal component analysis of the 
light curves can yield a less biased infor mation about the 
mass function of the lenses as shown bv lAlardI ll200lh . 
iHanl (J_996) has estimated the amount of expected mi- 
crolensing events with a cumulative event signal- to-noise 
ratio (S'/iV)min — 20 (for a further description see ICouldl . 
1996). The assumption he used was that the whole dark 
matter halo consists of MACHOs. The calculations are 
mainly depending on the luminosity function of the M31 
stars and on the spatial distribution of the lenses. From 
Han's paper we derive the following scaling relation for 
the event rate: 



r = 94ev/yr 
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/ ^scc \ ^ / ^sky + ^psf \ / ^ccd \ / 

\l" J [ 3.35 J UoarcminV ^' 

For the typical parameters of our 1999/2000 campaign 
this results in a predicted event rate F = 27ev/yr, with 
the length of the campaign icamp = 260 d, the average 
time span between observations tcyc = 2 d, the exposure 
time per night icxp = 0.5 h, the photon detection rate 
e = 43phot/s, the median seeing ^scc — 1.5arcsec, the 
area of the image Qccd ~ TOarcmin^, the inaccuracy of 
the PSF matching rj^^f = 0, and the factor in the noise 
due to sky photons rj'^j^y — 2.6. 

This calculation can only serve as a crude estimate of the 
event rate. In our subsequent papers we will present cal- 
culations more appropriate for our data set. 

3. The WeCAPP Project 

3.1. Telescopes and Instruments 

The Wendelstein 0.8 m telescope has a focal length / of 
9.9 m, which results in an aperture ratio f/D — 12.4. 
Starting in September 1997 we used a TEK CCD with 
1024 X 1024 pixels of 24 /im corresponding to 0.5arcsec 
on the sky. With this CCD chip we were able to cover 
8.3 X 8.3arcmin^ of the bulge of M31. To increase the 
time sampling of our observations we started to use the 
Calar Alto 1.23 m telescope (/ = 9.8 m, f/D = 8.0) in 
1999. The observations were partly carried out in service 
mode. Six different CCD chips were used. Three of these 
CCDs cover a field of 17.2 x 17.2arcmin^ and were used 
to survey the whole bulge for lensing events. A detailed 
overview of the properties of each CCD camera used for 
WeCAPP is given in Table [H 

Most of the sources for possible lensing events in the bulge 
of M31 are luminous red stars i.e. giants and supergiants. 
Consequently the filters used in our project should be sen- 
sitive especially to these kind of stars. We chose there- 
fore R and I filters for our survey. At Wendelstein we 
used the R2 (A ~ 650 nm, AA ~ 150 nm) and Johnson I 
(A ~ 850 nm, AA ~ 150 nm) wavebands. To be as consis- 
tent as possible with the data obtained at Wendelstein, the 
Calar Alto observations were carried out with the equiva- 
lents filters, R2 (A ~ 640 nm, AA ~ 150 nm) and Johnson 
I (A ~ 850 nm, AA ~ 150 nm). Since June 2000 we are 
using the newly installed filters Johnson R (A ~ 640 nm, 
AA ~ 160 nm) and Johnson I (A ~ 850 nm, AA ~ 150 nm) 
at Calar Alto. 

Despite of the combination of different telescopes, CCDs, 
and slightly different filter systems we observed no sys- 
tematic effects in the light curves depending on these pa- 
rameters. 

3.2. Observing Strategy 

To follow the suggest ion of lTomanev fc: CrottsI l)l996() and 
iHan fc Could* ('1996') we chose the field with the maxi- 
mal lensing probability, pointing to the far side of the 
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Fig. 1. M31 frame observed on 2000 June 26th with the 
Calar Alto 1.23 m telescope: 17.2 x 17.2 arcmin^ (3.85 x 
3.85 kpc^). A quarter of this field, marked by the solid 
lines, was observed during the first three campaigns from 
September 1997 until March 2000. 



M31 disk. The main fraction of the field is covered by the 
bulge of M 31 with the nucleus of M 31 located at one cor- 
ner of the field (Fig.^. As gravitational lensing is achro- 
matic, the amplification of the source is the same in differ- 
ent waveband s . However, as shown in several pape rs (e.g. 
'Va,]ls-Cab aiidl ll 994 Iwitl ll 9951 iHan et a,].L l2000h blend- 
ing on the one hand and differential amplification of an 
extended source on the other hand can lead to a chro- 
matic, but still symmetric, lensing light curve. Under cer- 
tain circumstances chromatic light curves permit to con- 
strai n the physical properties of the source-le ns system 
(e.g. ICould fc Welcil Il996l: iHan & Parkl l200lh . Variable 
stars will generally change colour in a different way. Our 
observation cycle therefore comprises 5 images in the R 
band and 3 images in the I band lasting about 45 min 
including readout time. Stacking these images with an av- 
erage exposure time of 150 sec in R and 200 sec in I results 
in a magnitude limit between (20.8 - 22.1) mag in R and 
(19.1 - 20.4) mag in I for a point source on the back- 
ground of M31 and a signal-to-noise ratio (S/N) — 10 in 
over 95 % of the frame. The background of M 31 typically 
has a surface brightness between (18.7 - 21.2) mag/arcsec^ 
in R and (16.8 - 19.3) mag/arcsec^ in I. The cycles were 
repeated as often as possible during one night, usually at 
least twice. As we had to avoid saturation of stars in the 
observed field we made exposure times dependent of the 
actual seeing, whereas exposure times in the I bands where 
generally longer. 
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Site 


Campaign 


CCD 


Size 


[arcsec /px] 


Field [arcmin^] 


# of R frames 


# of I frames 


We 


1997/1998 


TEK#1 


IK X IK 


0.49 


8.3 X 8.3 


276 


123 


We 


1998/1999 


TEK#1 


IK X IK 


0.49 


8.3 X 8.3 


454 


210 


We 


1999/2000 


TEK#1 


IK X IK 


0.49 


8.3 X 8.3 


835 


358 


CA 


1999/2000 


TEK7c_12 


IK X IK 


0.50 


8.6 X 8.6 


266 


136 


CA 


1999/2000 


TEK13C.15 


IK X IK 


0.50 


8.6 X 8.6 


62 


33 


CA 


1999/2000 


SITe2b_ll 


2K X 2K 


0.50 


17.2 X 17.2 


14 


7 


CA 


1999/2000 


SITc2b_17 


2K X 2K 


0.50 


17.2 X 17.2 


448 


249 


CA 


1999/2000 


SITel8b_ll 


2K X 2K 


0.50 


17.2 X 17.2 


165 


91 


CA 


1999/2000 


LORlli_12 


2K X 2K 


0.31 


10.75 X 10.75 


219 


92 



Table 1. Properties of all CCD cameras used during WeCAPP at Wendelstein (We) and Calar Alto (CA) 
Observatories, respectively. All CCDs have a pixel size of 24 /um, except the Loral which has 15 /xm pixels. 



Filter 


We 1997/1998 


We 1998/1999 


We 1999/2000 


CA 1999/2000 


R 


2.76 


1.45 


1.40 


1.49 


I 


2.59 


1.44 


1.32 


1.44 



Table 2. Median values of the FWHM of the PSF, given in arcsec, for the images taken during WeCAPP in the R 
and I band at Wendelstein and Calar Alto Observatories. 



120 1 




PSF PSF 

Fig. 2. Histograms of the Full Width Half Maximum (FWHM) of the point spread function (PSF) of the frames taken 
at Wendelstein Observatory during the 1997/1998 campaign (left panel) and the 1998/1999 campaign (right panel). 
Frames in the R band are marked by a solid line, frames in the I band by a dashed line. The lower limit of the PSF is 
restricted by a pixel size of 0.5 arcsec. 
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Fig. 3. Histograms of the FWHM of the frames taken during the 1999/2000 campaign at Wendelstein (left panel) and 

Calar Alto Observatory (right panel). Frames in the R band arc marked by a solid line, frames in the I band by a 
dashed line. Note that the pixel sizes of the CCD cameras used correspond to 0.5 arcsec on the sky. 
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Julian date - 2451000 

Fig. 4. Illustration of PSF vs. time coverage during three years of WeCAPP. Shaded regions mark the periods of time 
when M31 was not observable. 



3.3. The Data 

We began our observations at Wendelstein with a test 
period in September 1997, observing on 35 nights until 
March 1998. The second observational period lasted from 
1998 October 22nd until 1999 March 24th . During the 
first Calar Alto campaign we received two hours of service 
observations on 87 of 196 allocated nights (1999 June 27th 
- 2000 March 3rd). From November 1st until November 
14th we were able to observe during the whole night. In 
parallel we continued our observations at Wendelstein on 
221 nights, of which 65 were clear. In this way we achieved 
an overall time coverage of 132 nights (52.6%). During 
the three years of WeCAPP we collected at Wendelstein 
Observatory a total of 1565 images in the R band and 691 
images in the I band. The observations for over more than 
one year at Calar Alto resulted in 1174 frames in the R 
and 608 frames in the I band. 

During the 1997/1998 test campaign conditions at the 
Wendelstein telescope were improved significantly. A new- 
ly installed air conditioning system reduced dome seeing 
to a low level. Further improvements like fans just above 
the main mirror finally lead to a leap in the image quality 
obtained with the telescope. Figure El which presents the 
PSF statistics of Wendelstein images from the 1997/1998 
and 1998/1999 campaigns respectively illustrates this fact. 
In general Wendelstein shows a marginally better PSF dis- 
tribution than Calar Alto (see Figs. El and OJ. Table (21 
shows the PSF median values for the images taken during 
WeCAPP at both sites. 

Figure^ shows the time sampling we reached so far with 
WeCAPP. Because of time loss during the upgrades of 
the telescope, time coverage of the 1997/1998 campaign is 
only fragmentary. About the same applies to the follow- 
ing campaign, this time due to a camera shutdown and 
another time consuming project. Finally time coverage of 
the first joint campaign of Wendelstein and Calar Alto is 
good, last but not least due to the often opposite weather 
situation in Spain and Germany. 



iGossl fc RiffesCTi l|200lh and combines all reduction steps 
from de-biasing of the images until the final measurements 
of the light curves in one software package, including full 
error propagation from the first reduction step to the last: 

1. standard CCD reduction including de-biasing, flat- 
fielding and filtering of cosmics 

2. position alignment using a 16 parameter interpolation 

3. stacking of frames 

4. photometric alignment 

5. PSF matching using OIS ( Optimal Image Subtra c- 
tion), a method proposed bv I Alard fc LuDtonI l(l998l) 

6. generation of difference images 

7. detection of variable sources 

8. photometry of the variable sources 

4.2. Standard CCD Reduction 

Pre-reduction of the raw frames is performed in a stan- 
dard way: After de-biasing of the frames, saturated and 
bad pixels are marked. We use the 3 a clipped median of 
a stack of at least 5 twilight flatfields for the flat fielding 
procedure. 

Cosmics are effectively detected by fitting a Gaussian PSF 
to all local maxima in the frame. All PSFs with a FWHM 
lower than 0.7 arcsec and an amplitude of 8 times the 
background noise are removed. This provides a very reli- 
able identification and cleaning of cosmics. 

4.3. Position Alignment 

After determining the coordinates of the reference objects 
by a PSF fit, we calculate a linear coordinate transforma- 
tion to project an image onto the position of the reference 
frame. 

The flux interpolation for non-integer coordinate shifts is 
calculated from a 16-parameter, 3"^-order polynomial in- 
terpolation using 16 pixel base points. 



4. Data Reduction 

4.1. Reduction Pipeline 

In the last three years we developed an image reduction 
pipeline being able to cope with a massive imaging cam- 
paign. This reduction pipeline is described in detail in 



4.4. Stacking Frames 

To avoid saturation of Galactic foreground stars and the 
nucleus of M31 in our field, exposure times were limited to 
a few hundred seconds. Therefore one has to add several 
frames taken in one cycle to obtain an acceptable signal- 
to-noise ratio (S/N). Usually we stacked 5 frames in the R 
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band and 3 frames in the I according to two criteria, com- 
parable PSF and comparable sky. Frames with very high 
background levels or very broad PSFs were not added if 
they reduced the detectability of faint variable sources. 
Consequently the number of images to be stacked was not 
fixed, coaddition of frames was performed in a way to 
get a maximum (S/N) ratio for faint point sources in the 
stacked frame. 

Before OIS is finally carried out we align all frames pho- 
tometrically. This ensures that all light curves are photo- 
metrically calibrated to a standard flux. 

4.5. PSF Matching 

In order to extract light curves of variable sources from 
the data we use a m ethod called Difference Image Analysis 
(DIA), pro posed bvlCiarduUo et all lIlQQOl) and first imple- 
mented bv lXomanev &: CrottTijTQQGl) in a lensing study. 
The idea of DIA is to subtract two positionally and pho- 
tometrically aligned frames which are identical except for 
variable sources. The resulting difference image should 
than be a flat noise frame, in which only the variable point 
sources are visible. 

The crucial point of this technique apart from position reg- 
istration is the requirement of a perfect matching of the 
point spread functions (PSFs) between the two frames. 
The PSF of the reference frame r is convolved with a ker- 
nel k to match the broader PSF of an image i, 



i(x, y) ~ r{u, v) k{u, v) + bg{x, y) = f{x, y) 



where bg accounts for the background level and f is the 
transformed reference frame. 

In order to obtain a n optimal kernel k we implemented 
OIS as proposed bv lAlard fc LuptonI l|l998l) . This least- 
squares fitting method determines k by decomposing it 
into a set of basis functions. We use a combination of 
three Gaussians with different widths cr multiplied with 
polynomials up to 6*'' order. This leads to the following 
49 parameter decomposition of k{u, v): 



CTi = 1 : e '"t (ai + 

,,2 I ,,2 



a28^''^) 



(T2 = 3 



'I" (a29 + ... + asgu'' + ... + 043^"') 



(6) 



(Ta = 9 : e ^"3 (044 + . . . + qaiu^ + a4suv + 04911^) 
Additionally 3 parameters are used to fit the background 



bg{x, y) = + a^ix + 052 y . 



(7) 



To cope with the problem of a PSF varying over the area 
of the CCD we divide the images in sub-areas of 141 x 
141 pixels each. In all sub-areas a locally valid convolu- 
tion kernel is calculated. As we have chosen the kernel to 
have 21 X 21 pixels we therefore effectivly use 161 x 161 
to derive k{u, v). 

Differentia l refraction causes a star' s PSF to depend on 
its colour i Tomanev fc Crott3 . Il996 l. chap. 4.4). However 
these second order effects are negligible for our data set 



and do not lead to residuals in the difference images. 
As we are performing DIA we have to choose a reference 
frame r which will be subtracted from all other coadded 
frames i and which determines the baseline of the light 
curve. OIS shows best results for a small PSF and a high 
{S/N) reference frame. Therefore, the best stacked im- 
ages were coadded once more. Our actual R band refer- 
ence frame comprises 20 images taken at 2 different nights 
resulting in a total exposure time of 2400 seconds and a 
PSF of a FWHM of 1.05 arsec (= 2 pixels for aU CCDs ex- 
cept one at Wendelstein and Calar Alto). For the I band 
we coadded 18 frames (i.e. 6 stacked frames taken at 3 
different nights) which results in a total exposure time 
of 3160 seconds and a FWHM of 1.15 arcsec. As we are 
continuing collecting data the process of constructing the 
ultimate reference frame has not finished yet. Each night 
of high quality data collected at one of the two sites will 
improve the reference frame further. Figure|Slshows a typi- 
cal difference image obtained by using our implementation 
of OIS. 

4.6. Source Detection 

To detect source s in the differe nce images we fit a rotated 
Moffat function (iMoffad Il96l to all local maxima in the 
binned frame. Real sources are filtered by rejecting sources 
with an amplitude less than 5 times the background noise. 



(5) 4. 7. Photometry of the variable sources 



Photometry of the detected sources is performed by a pro- 
file fitting technique. To obtain information about t he PSF 
of an y particular frame we apply a Moffat fit (fMoffatl . 

to several reference stars in the CCD field. Having 
determined the shape of the PSF, we perform Moffat fits 
on the positions of the variable sources as returned from 
the detection algorithm. In these final fits the amplitude 
is the only free parameter. To determine the fiux of the 
source we finally integrate the count rates over the area of 
the (now fully known) analytical function of the PSF. This 
minimizes the contamination from neighbouring sources. 

4.8. Calibration of the Light Curves 

As the coadded images are normalized to the reference 
frame it is not necessary to calibrate each image sepa- 
rately. Only the reference frame is calibrated once. 
To calculate magnitudes tor/ in our R' band, which corre- 
sponds to R2, we determined the instrumental zcropoint 
ZP^i' and the extinction coefficient k: 



Tor/ 



-2.5 \og[ADU/t\ - K • AM + ZP^ 



(8) 



where t is the exposure time and AM is the airmass. 
Apert ure photometry with 7 different Landolt standard 
stars llLandold. Il992l) observed at different airmasses 
was performed for a photometric night at Calar Alto 
Observatory. With these stars the extinction k for the 
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Fig. 5. Difference images of a part of tlic M31 bulge (3.25 x 2.82 arcmin^), centered at 6.6 arcmin distance from the 
nucleus. Left panel: difference frame without OIS. Because of the relatively large residuals from stars no identification 
of faint variable stars is possible. Right panel: difference frame with OIS. Bright and dark spots are variable sources. 
The two boxes represent stars not subtracted in order to provide information on the PSF. 



night was calculated to kr' = 0.073 ±0.005. To determine 
the zeropoint for the R' band we used an AOV-star, Feige 
16, with the colours (B -V)^ -0.012, {U - B) = 0.009, 
(y - i?) = -0.003, (R- I) ^ 0.002, and a visual magni- 
tude oi V — 12.406 mag. The zeropoint was determined 
according to Eq. |Hlto ZPr' = 23.05 ± 0.02 mag and used 
to calculate the magnitudes for the reference frame. This 
zeropoint is not valid for Wendelstein. 
In the following, we only give fluxes for the sources in our 
filter system, because the intrinsic magnitudes and colours 
of our unresolved sources cannot be determined with suf- 
ficient accuracy. 

We show the light curves in flux differences according to 



^/r' — /R'.Vcga 



AADUn' 
t 



10' 



0.4 k am 



10 



-0.4 ZP„ 



(9) 



with /R'.vcga = 3124 Jy from an integration over the CCD- 
filter system. 

The same transformations were done for our I' band, cor- 
responding to Johnson I (Calar Alto), with Kp = 0.025 ± 
0.005, ZPv = 21.82 ± 0.03 mag and /I'^voga = 2299 Jy. 
The colour terms between all filter sets we used in our ob- 
servation are negligible. The transformation to the stan- 
dard Kron-Cousins filter system is: R = R' + ZPy<' + 
OM{R -I),I = r + ZPv + 0.38(i? - /). 

5. Results 

We present a small sample of light curves to show the effi- 
ciency of the method. All light curves were observed over 
more than three years from 1997 until 2000. Time spans 
when M 31 was not observable are marked by shaded re- 
gions. Because of bad dome seeing conditions and an in- 
appropriate autoguiding system errors were largest during 
the first Wendelstein campaign 1997/98. During the sec- 
ond period 1998/99 we were able to decrease the FWHM 
of the PSF by a factor of two, thus the photometric scatter 




-2 2 
ADU/pixel 

Fig. 6. Histogram of the pixel values of a simulated dif- 
ference image dxy = ixy — fuv ^ kuv {r is the reference 
frame, i is a second frame and k is the convolution kernel). 
The pixel values are divided by the expected rms errors, 
as derived from error propagation. The solid curve shows 
a Gaussian with a ~ 1. We calculated the reduced chi- 
square Xdof 19 different simulated images in the range 
between -3 and 3. The median Xdof l-^' which means 
that expected and measured errors match almost perfectly 
and that the residuals of the OIS are at the photon noise 
level. 



is also clearly smaller. During the third period 1999/2000 
we observed simultaneously at Calar Alto and Wendelstein 
and got data points for 53 % of the visibility of M31. 
The OIS method can be applied for very crowded fields 
like M31 and gives re sidual errors at the ph oton noise 
level (Fig. EJ, see also lOossl fc Riffeseii lj200l|) V A good 
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Fig. 8. Light curve of the (5-Cephei star of Fig. din the R' band, convolved with its period of P = 15.76 ± 0.01 d. 
Plotted without (left panel) and with (centre panel) la error bars, which represent fully propagated errors through 
all reduction steps. Right panel: Binned R' light curve of this star. 



estimate for the average noise present in the area flpsF 
of a PSF is iV = 0.1 x 10~^Jy. The hght curves of vari- 
able stars presented in the Figs. 171 through [TCI indicate a 
typical scatter which is in good agreement with the above 
estimate. This means that a red clump giant with a bright- 
ness of Ml = jGrillmalr et al, . 1996. Fig. 7) and a colour 
of {R — I) — 0.5 |Leieune et aL . ll998t) has to be amplified 
by a factor of 10 to be detected with a peak signal-to-noise 
ratio of {S/N) = 3 in our survey. The brightest RGB stars 
with a Ml = —3.5 and a colour of (i? — /) = 1 need an 
amplification of 1.6 only. 

Up to now we detected over 5000 variable sources in a 
8x8 arcmin^ field. A preliminary analysis of the light 
curves shows that we have found the whole range of vari- 
able stars including novae and other types of eruptive vari- 
ables, Cepheids, semi-regular, Mira-type and other long- 
periodic variables. 

In Fig. [3 we present one of the (5-Cephei variable stars in 
the R' and 1' bands, Fig. shows the R' light curve of 
this star convolved with its period, which was determined 
to 15.76 ± 0.01 days. Figure El presents the light cur ve of 
a nova previously pubhshed bv lModiaz fc Lil l)l999(l . It's 
the brightest variable source detected in our M31-field. 
Figure IIUI is an example for an eruptive variable star, 
which could be mistaken as a microlensing event, if the 
time coverage were insufficient. Figures ITTI to [T51 display 
light curves of variable stars, which were classified as long- 
periodic in a preliminary analysis. Finally we present the 
light curve of a RV Tauri star in Fig. ^1 



6. Conclusions 

We presented an overview of the Wendelstein Calar Alto 
Pixellensing Project (WeCAPP). We demonstrated that 
despite observing at different sites with different instru- 
ments all dat a can be used fo r optim al image subtrac- 
tion following lAlard fc LuptonI l)l998|) . This method can 
be applied for very crowded fields like M 31 and gives 
residual errors at the photon noise level. A red clump gi- 
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Fig. 7. Light curve of a (5-Cephei variable, upper panel: 
R' band, lower panel: 1' band. 
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Fig. 9. Light curve of a nova, representing the brightest 
variable source detected in our M 31-field. This nova was 

Lil l|l999^ . Upper panel: 



previously published bv lMoaiaz 
R'-Band, lower panel: I'-Band 
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Fig. 10. Light curve of an eruptive variable, which could 
be mistaken as a microlensing event, if the time coverage 
were insufBcient. Upper panel: R' band, lower panel: I' 
band. 
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Fig. 12. Light curve of a longperiodic variable. Upper 
panel: R' band, lower panel: I' band. 
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Fig. 11. Light curve of a longperiodic variable. Upper 
panel: R' band, lower panel: I' band. Note, that insuffi- 
cient time coverage could result in a false identification of 
this variable as a microlensing event. 

ant of M\ = 0, which is amplified by a factor of 10 by 
a microlensing event, can be detected with our data. Wc 
showed how the data are reduced and how light curves are 
extracted. For illustration we presented a small sample of 
light curves. In future publications wc will present a full 
catalogue of variable sources which we found in our M 31 
field, including potential MACHO light curves. 
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Fig. 13. Light curve of a longperiodic variable. Upper 
panel: R' band, lower panel: I' band. 
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Fig. 14. Light curve of a longperiodic variable star with a 
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lower panel: I' band. 
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Fig. 15. Light curve of a longperiodic variable. Upper 
panel: R' band, lower panel: I' band. 
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Fig. 16. Light curve of a RV Tauri star in the R' (upper 
panel) and I' (lower panel) bands. Due to the optimal time 
coverage, the typical double- wave shape with alternating 
deep and shallow maxima of the light curves of this class 
of variable stars is uncovered. 
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